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, , Abstract 

, Under assumption of singular behavior of as{q^) at ~ and of large behavior, corre- 

' spending to the perturbation theory up to four loops, a procedure is considered of matching the 

/3-function at a boundary of perturbative and non-perturbative regions. The contribution of the 
non-perturbative region to the gluon condensate is calculated with varying normalization condition 
(N . {ml) = 0.29, 0.30, 0.36 for two different ways of definition of the non-perturbative invariant 

charge in the infrared region. We obtain quite consistent results for values of the gluon condensate, 
nonperturbative region scale qo, and the string tension a. 
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It is well-known, that the /3-function in the perturbative QCD is of the form 



For Uf — 3 we have values of coefHcients &o = 9, 6i = 32, 62 — 1287.67, 63 ~ 12090.38 (coefficients 60, 
bi do not depend on renormalization scheme while values 62, 63 correspond to a choice of MS'-scheme). 
Expressions obtained by solution of Gell-Mann - Low equation 



Oh! 

^: = (2) 



with the use of (1), are widely used for sufficiently large momenta transfer, however they can not be 
^ , applied in the infrared region. For a behavior of the invariant charge (the running coupling constant) 

' as{q^) at 5^ — > a number of variants are considered (see, e.g., [1]). In particular singular infrared 

asymptotic at = 



52 m2 
47r (j2 



is possible (see, e.g., review [2] and more recent papers [3]). Such behavior naturally corresponds to a 
linear confining quark-antiquark static potential at that g^M"^ = Gttit, where a is the string tension. 
Results of some works on the lattice study of the three-gluon vertex [4] demonstrate a necessity of taking 
into account of non-perturbative contributions to the running coupling constant being of the form of (3). 
In the framework of the continuous QFT arguments in favor of behavior (3) are also presented in recent 
paper [5]. 

Asymptotic behavior (3) occurs provided —t —h for h 00. Let us consider a possibility of 
behavior (3) to be valid and assume the following form of the infrared /3-function 

f3{h) = -h + z, h>ho, (4) 
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where z is a constant and /iq defines the boundary between perturbative and non-perturbative regions. 
For h < ho we shah use /3-function (1) with a finite number of loops taken into account. Our recipe for 
construction of /3-function for all /i > consists in a smooth matching of expressions (4) and (1) at point 
h = ho in approximations of the perturbation theory up to four loops. The demand of the /3-function 
and its derivative to be continuous uniquely fix free parameters z and ho of the "global" /3-function (the 
matched one) . 

For an illustration let us consider the most simple one-loop case. Conditions of matching give two 
equations 

-bohl = -ho + z, 
-2boho = -1. (5) 

The solution of set (5) reads 

We shall normalize perturbative solution 

/ 2n 47r g2 

by value 47r/io, that gives 

2^0 = ql/^QCD = e^, (8) 

where e = 2.71828... . Imposing on as(g^) the natural condition to be continuous at = Qq, we may 
normalize non-perturbative solution of equation (2) 



a,{q^) = 47r{^ + z], < ql (9) 



by 47r/io as well. As a result we obtain 



C = ql{ho- z), co = C/h?QCD = xo{ho-z) = -^^. (10) 

For final fixation of the solution for all > we need to define Aqcd by normalizing the solution, say, 
at point = m^, where rrir = 1.77703 GeV is the mass of the r-lepton [6]. 

Let us turn to calculation of the gluon condensate. Its value is defined by the non-perturbative part 
of as- We have (see, e.g., the third of refs. [3]) 

oo oo 

K =< as/n : Gl^G^, :>= ^ j dq^ ?'a„pt(g') = ^ J ^q^ q^Mq^) " apert(g')). (11) 



In our approach the non-perturbative contribution is present for < only. For the beginning we 
define the perturbative part at this region basing on the assumption of freezing of at small q^ (see, 
e.g., [7]). That is we assume 

Q!pert(9^) = Ois{ql) =47r/io, q^<ql. (12) 
Using expressions (9), (12), we have 

il 

I^='^2jdq'q'{j,+z-ho)^ f,qt (| + ^(^ " ^^o) 
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{ho - z)xlK%cD- (13) 
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Expression (13) is valid for each of 1 — 4-loop approximations of the perturbative and ratio K/ Aq^jjj 
does not depend on a normalization of a^. In the one-loop case we have from (6), (8), (13) with bo = 9 

c4 



(14) 



In the one-loop case normalization condition as(m^) = 0.32 gives Aqcd = 0.201 GeV, and expression (14) 
gives K = (0.197 GeV)^, which is far from the conventional value of the gluon condensate (0.33 ± 0.01 

GeV)4 [8]. 

Let us consider multi-loop cases. Solution h{q^) of equation (2) for L — ln(q^/A-^) oo reads as 
follows 



1 
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4bi 
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In i-lnL - 1 + 
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(15) 



Keeping in the expression terms with powers of logarithms in denominators up to the first, the second, 
the third and the fourth, we fix the 1 — 4-loop approximations of the perturbation theory for running 
coupling constant. It may be written in the form 



a,(g2) =47r/i(g2) 



47r 
6^ 



a{x), 



I \ 1 ,ln(lna:) 

aix) = b \ + b^ 

Inx In^x 



ln^(lna;) In (Inx) k 



In^x 



ln^(lna:) 51n^(lnx) ln(lna;) 



In"* X 



2 In^x 



In* a; 



Here x = J , and coefficient are defined as follows 



+ 



In^a; 



In^x 



(16) 



K = 



26i 

"o 

-1 + 



6062 
86f 



166? ■ 



(17) 



Coefficients 6, k, k depend on n/.With n/ = 3 we have 6 ~ 0.7901, k ~ 0.4147, k ~ —1.3679. In the case 
of the two-loop approximation for perturbative as we have the following set of equations 



The set for three loops reads 



60/10 + 261/10 = ho - z, 
260/10 + 661/1^ = 1, 



ln(lna;o) 
b—^ = boho, 



1 



lna;o 



£0 

Xo 



In^ Xo 
+ z = ho. 



(18) 



6o/i^ + 26i/i^ + y/i^ = /io--2, 
260/10 + 6bihl + 2b2hl = 1, 
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In a;o 



ln(ln .To) 
In^ Xq 



In^(lnxo) ln(lna;o) 



In^ .To 
Co 



In xq In a;o 



boho, 



(19) 



Xo 



The set for four loops reads 



+ z = ho- 



boh'o + 2bih^ + jh^ + b^h^ = ho-z, 
2boho + 6bihl + 262/10 + S^s/io = 1> 



In Xo 



ln^(lna;o) ln(lna;o) 



ln(ln To) 
In^ Xo 



In^(lnTo) 5 In (lna;o) 
In^ Xo 2 In^ xq 

Co 
Xo 



In^ Xo 
+ {3k + 1) 



In'^xo 



In^ Xo 

In(lnTo) ^ K 
In^ Xo In^ Xo 



= boho, 



(20) 



+ z = ho- 



From sets of equations (18) - (20) we find values of ho, z, Xo and o- They are presented in Table 1. 
The value of K^^'^/Aqcd calculated with the aid of expression (13) is also presented there. For the 
parameter of string tension a by using expressions (9), (10) we have relation cr/h^^cD = Sttcq/S. Taking 
into account the existing data [6, 9, 10] we fix the momentum dependence of solutions for a number of 



Table 1: The dimensionless parameters ho, z, xo = Qo/^qcd' Co = C/K^qjj, K^/^ / Kqcd on the number 
of loops, n/ = 3 





1-loop 


2-loop 


3-loop 


4-loop 


ho 


0.0556 


0.0392 


0.0356 


0.0337 


z 


0.0278 


0.0215 


0.0203 


0.0197 


Xo 


7.3891 


7.7763 


10.2622 


12.4305 


Co 


0.2053 


0.1374 


0.1572 


0.1741 




0.9799 


0.8977 


0.9952 


1.0709 



values of the running coupling constant at r-lepton mass scale rrir with the effective number of flavors 

Uf being equal to 3. The corresponding to these normalization conditions values of Aqcd are presented 
in Table 2, values of boundary momentum qo = ^/xqAqcd are presented in Table 3, values of the gluon 
condensate K^/^ are presented in Table 4 and the string tension parameter a is given in Table 5. 



Table 2: Values of the parameter Aqcd (GeV) on loop numbers and normalization conditions. Normal- 
ization conditions: a^m^) = 0.29. 0.30 0.36 with m.r = 1.77703 GeV. n/ = 3 





1-loop 


2-loop 


3-loop 


4-loop 


0.29 


0.1600 


0.3168 


0.2873 


0.2837 


0.30 


0.1734 


0.3370 


0.3069 


0.3026 


0.31 


0.1869 


0.3568 


0.3263 


0.3212 


0.32 


0.2005 


0.3762 


0.3454 


0.3394 


0.33 


0.2143 


0.3951 


0.3642 


0.3573 


0.34 


0.2280 


0.4136 


0.3827 


0.3749 


0.35 


0.2418 


0.4315 


0.4007 


0.3920 


0.36 


0.2556 


0.4490 


0.4184 


0.4087 
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Table 3: Values of the parameter qo (GeV) on loop numbers and normalization conditions. Normalization 
conditions are the same as in Table 2 



S V T / 


1-loop 


2-loop 


3-loop 


4-loop 


0.29 


0.4350 


0.8833 


0.9203 


1.0002 


0.30 


0.4713 


0.9397 


0.9831 


1.0667 


0.31 


0.5081 


0.9949 


1.0452 


1.1323 


0.32 


0.5451 


1.0490 


1.1065 


1.1967 


0.33 


0.5824 


1.1018 


1.1667 


1.2598 


0.34 


0.6198 


1.1533 


1.2258 


1.3216 


0.35 


0.6572 


1.2034 


1.2838 


1.3820 


0.36 


0.6947 


1.2522 


1.3405 


1.4409 



Table 4: Values of the gluon condensate K^^'^ (GeV) on loop numbers and normalization conditions. 
Normalization conditions arc the same as in Table 2 





1-loop 


2-loop 


3-loop 


4-loop 


0.29 


0.1568 


0.2844 


0.2859 


0.3038 


0.30 


0.1699 


0.3025 


0.3054 


0.3240 


0.31 


0.1832 


0.3203 


0.3247 


0.3439 


0.32 


0.1965 


0.3377 


0.3437 


0.3635 


0.33 


0.2099 


0.3547 


0.3624 


0.3827 


0.34 


0.2234 


0.3713 


0.3808 


0.4014 


0.35 


0.2369 


0.3874 


0.3988 


0.4198 


0.3(i 


0.2504 


0.4031 


0.41()4 


0.4377 



Table 5: String tension parameter a = yja (GeV) on loop numbers and normalization conditions. Nor- 
malization conditions are the same as in Table 2 



as(m^) 


1-loop 


2-loop 


3-loop 


4-loop 


0.29 


0.2098 


0.3398 


0.3297 


0.3426 


0.30 


0.2274 


0.3615 


0.3522 


0.3654 


0.31 


0.2451 


0.3827 


0.3745 


0.3878 


0.32 


0.2630 


0.4035 


0.3964 


0.4099 


0.33 


0.2809 


0.4239 


0.4180 


0.4315 


0.34 


0.2990 


0.4437 


0.4392 


0.4527 


0.35 


0.3170 


0.4()29 


0.4599 


0.4733 


0.36 


0.3351 


0.4817 


0.4802 


0.4935 



Let us consider another variant of definition of the perturbative ols behavior for < q^. Namely 
instead of freezing (12) we shall assume "forced" analytic behavior of as in this region, 

apert(g^) = Q;an(9^), 9^ < Qq- (21) 

The main ideas of the analytic approach in quantum field theory, which allows one to overcome difficulties, 
connected with nonphysical singularities in perturbative expressions, were proposed in works [11, 12]. The 
analytic approach is successfully applied to QCD [13]. The forced analytic running coupling constant is 
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defined by the spectral representation 



aan(y) 



(22) 



where spectral density p{a) = 9aan(— f — ^0) = S>a(— cr — iO). For the perturbative solutions a{x) of the 
form (16) the two- loop analytic running coupling constant and its non-pcrturbativc part were studied in 
ref. [14], the three-loop case and the four-loop case where studied in refs. [15] and [16], respectively. Let 
us write the spectral density up to the four-loop case. 



p(^)(a)=pW(a) 



i2 + 7r2 ' 



[27rtFi{t) - {t-" - n^) F2{t)] , 



[tt (3*2 - TT^) {F^{t) - Fi{t)) - 2t - 377^) Fi{t)F2{t) 



Here t = ln((j). 



(i2 + 7r2)3 

-TT (3*2 - TT^) Fi{t) +t [t^ - 377^) F2{t) + TTK (3*^ - TT^)] , 

pW(a) = p('\a) ^2^, 6.2,2 ^ ^4) (^3(,) _ 3i^f (t)i^.(t)) 

+477* (*2 - TT^) (Ff (*) - 3^^l(*)F|(t)) - lO^t (*2 - (^2(,) _ ^2(,)) 

+5 (*■* - 677^*2 + tt'') Fi{t)F2{t) + 47r (1 + 3k) * ~ tt^) Fi(*) 
- (1 + 3k) (*" - 67r2*2 + tt") f2(*) + 47rK* (*2 - tt^)] . 



Fi(*) = i ln(*2 + 7r2), i^2(*) = arccos ^^p=f • 



(23) 
(24) 

(25) 



(26) 
(27) 

(28) 



Solving the equation^ 

aan(?/) =boho, 

we find values yo for Oan, being defined by formulas (22) - (27) with the use of values ho obtained above. 
Further wc find dimensionless quantity ^ = {^an/ ^qcdY = ^o/vo- Values of yo and ^ in dependence 
on the number of loops are presented in Table 6. In Tabic 7 values of Aan = Qo/y/yo are presented in 
dependence on the number of loops and on normalization conditions at = m^. 

Let us turn to the gluon condensate. In the considered method of definition of the perturbative part 
of as in the non-perturbative region we have 



^an = J dq^ q^{Ois{q^) - Q:an(g2)) = ^ j dq^ + Z - ^aan ' 



(29) 



Taking into account the spectral representation (22) and performing integration in Eq. (29), we obtain 

oo 

12 At 



^an = ^ {Cql + ^qi) - ^ / da p{a) [yo - aln (l + . 



(30) 



Substitution a = cxp(*) with the use of expressions (8) and (10) leads to the following representation of 
the gluon condensate 



^an — ^%CD 



2 ^ 

^ (^0 - I) - 1^ y" dt Pit) {yoe* - e^* In (l + yoe"') } 



(31) 



^ While solving this equation it is convenient to use the method of ref. [15], in which aaS*(2/) is represented as a series in 

inverse powers of y 
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The expression in the square brackets of (31) does not depend on values of ag{m.'^), vahics of the ratio 
kH'^/Aqcd, which are obtained by numerical integration in formula (31), are given in Table 6. In Table 7 
and Table 8 values of the parameter Aan and of the gluon condensate kH"^ are presented, respectively. 

Table 6: Dimensionless parameters yo, kH'^ /Aqcd on loop numbers, n/ = 3 





1-loop 


2-loop 


3-loop 


4-loop 




1. ()()()() 


0.8299 


2.2091 


2.8710 




7.3893 


9.3702 


4.5225 


4.3297 




0.9373 


0.8494 


0.9353 


1.0025 



Table 7: Dependence of the parameter Aan (GeV) on loop numbers and normalization conditions. Nor- 
malization conditions are the same as in Table 2 





1-loop 


2-loop 


3-loop 


4-loop 


0.29 


0.43')() 


().9(i9() 


().()109 


0.5903 


0.30 


0.4714 


1.0315 


0.6526 


0.6296 


0.31 


0.5081 


1.0921 


0.6939 


0.6682 


0.32 


0.5451 


1.1515 


0.7345 


0.7063 


0.33 


0.5824 


1.2094 


0.7745 


0.7435 


0.34 


0.6198 


1.2659 


0.8138 


0.7800 


0.35 


0.6572 


1.3210 


0.8522 


0.8156 


0.36 


0.6947 


1.3745 


0.8899 


0.8504 



Table 8: Values of the gluon condensate Kan (GeV) in dependence on loop numbers and normalization 
conditions. Normalization conditions are the same as in Table 2 





1-loop 


2-loop 


3-loop 


4-loop 


0.29 


0.1500 


0.2691 


0.2687 


0.2844 


0.30 


0.1625 


0.2862 


0.2870 


0.3033 


0.31 


0.1752 


0.3031 


0.3052 


0.3219 


0.32 


0.1880 


0.3195 


0.3230 


0.3403 


0.33 


0.2008 


0.3356 


0.3406 


0.3582 


0.34 


0.2137 


0.3513 


0.3579 


0.3758 


0.35 


0.2266 


0.3666 


0.3748 


0.3929 


0.36 


0.2395 


0.3814 


0.3914 


0.4097 



For cases from one loop up to four loops the behavior of the running coupling constant as{q^) for all 
> is shown in Fig. 1. Here the behavior of the analytic coupling constant aa,n{q^) for q"^ < q^, which 
defines the perturbative part of asiq"^) in this region, is also shown. In Fig. 2 in addition to agiq"^) its 
non-perturbative part Q;npt('/'^) is also shown, which turns to zero for q^ > q^. Normalization condition 
for all curves in Fig. 1, Fig. 2 is as{m^) = 0.32. 

Starting from the well-known perturbative behavior of the /3-function for small values of the coupling 
constant h and the behavior for large values of the coupling constant, which corresponds to the linear 
confinement, we have constructed the /3-function for all h > 0. To control dependence of the results on 
the number of loops we simultaneously consider cases corresponding to perturbative /3-function for 1 — 4 
loops. While constructing the matched /3-function we assume it to be smooth at the matching point, that 
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cr 
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cr 
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Figure 1: Running coupling constant as{q^) for all q"^ > and analytic coupling constant aanC'Z^) for 
< (Jo (ths corresponding curves are the lower ones). Normalization conditions: as(m^) = 0.32, 



i(9o) = as(9o) 



leads to the invariant charge being smooth together with its derivative for all > 0. The normalization 
of the invariant charge, e.g., at fix it thoroughly. Then value of Q:s(m^) ~ 0.33 corresponds to string 
tension a ~ 0.42 GeV [17]. 

The obtained invariant charge is applied to study of a quite important physical quantity, the gluon 
condensate. In doing this we consider two variants of extracting of the non-pcrturbativc contributions 
from the overall expression for the invariant charge. The first variant assumes "freezing" of perturbative 
part of the charge in the non-perturbative region < Qq, while for the second one we choose the analytic 
behavior of the perturbative part in this region. As we sec from Tables 4, 8 the first variant leads to 
values of the gluon condensate being somewhat larger, than that for the second variant. Emphasize, that 
for as(m^) = 0.33 the gluon condensate for the second variant practically coincides with the conventional 
value [8], while for the first variant it turns to be K^^^ ~ 0.36 GeV, that is rather higher, than the 
conventional value. Note, that other important non-perturbative parameter go for the same normalization 
conditions also turns to be of reasonable magnitude, go — 1.17 GeV (Table 3). We may conclude. 
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Figure 2: Running coupling constant as^q^) and its non-pcrturbative part anpt(<7^) (the corresponding 
curves are the lower ones) with definition of apert(9^) by (21). Normalization condition: asim^) = 0.32 



that singular behavior of the running coupling constant (3) consistently describe these non-perturbative 
parameters. 

The authors express deep gratitude to V.A. Pctrov. V.E. Rochcv and D.V. Shirkov for valuable 
discussions. The work has been partially supported by RFBR under Grant No. 02-01-00601. 
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